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Abstract:- This paper studies a shock model for a repairable system with two-type failures by assuming that 
two kinds of shocks in a sequence of random shocks will make the system failed, one based on the inter-arrival 
time between two consecutive shocks less than a given positive value and the other based on the shock 
magnitude of single shock more than a given positive value. Further it is assumed that the system after repair is 
not 'as good as new', but the consecutive repair times of the system form a stochastic increasing a-series 
process. Under these assumptions, we determine an explicit expression for the average cost rate and an optimal 
placement policy N* based on the number of failure of the system is determined such that the long-run average 
cost per unit time is minimized. The explicit expression of long-run average cost per unit time is derived, and 
the corresponding optimal replacement policy can be determined analytically or numerically. Finally, a 
numerical example is given. 
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I. INTRODUCTION 

The shock model mainly used to study the external causes such as the environmental factor which may 
influence system reliability. For example, a precision instrument and meter installed in a power workshop may 
be affected by some random shocks due to the operation of other instruments, such as lathes or electrical 
machineries. Esary at al. [4] examined the theory of Poisson shock model. Later on, Barlow and Proschan [2] 
studied this problem in their monograph. Ross [14] developed a generalized Poisson shock model. Shantikumar 
and Sumita [16] extended Poisson shock to general shock, and studied a such shock model in which a system 
fails when the shock magnitude of single shock outstrips a given positive value. At the same time, Feldman [5], 
Zuckerman [24], Gottlieb [7] and Abdel-Hammed [1] determined respectively the optimal replacement policy 
for the different shock models. Li et al. [12] presented a new shock model named 5-shock model. Wang and 
Zhang [19] not only considered some reliability indices of the 5-shock model, but also studied geometric 
process repair model and its optimal replacement policy. 

For a deteriorating system, a more reasonable repair model is the geometric process repair model first 
proposed by Lam [8, 9]. Under this model, Lam considered two kinds of replacement policy, one based on the 
working age T of the system and the other based on the failure number N of the system. It assumed that the 
consecutive working times of the system form a decreasing geometric process while the consecutive repair times 
of the system form a increasing geometric process and the system after repair is not 'as good as new' state. 
Under these assumptions an explicit expressions for the long-run average cost per unit time under these two 
kinds of policy are respectively derived, and the corresponding optimal replacement policies T* and N* are 
found analytically or numerically. Because the geometric process is a special monotone process, Stadje and 
Zuckerman [17] introduced a general monotone process repair model to generalize Lam's work. Zhang [21] 
generalized Lam's work by a bivariate replacement policy (T, N) under which the system is replaced at the 
working age T or at the time of the N th failure, whichever occurs first. Later for improving the system reliability 
and its availability, Zhang [22] utilized the geometric process repair model to a two-component cold standby 
system with one repairman. It assumed that each component after repair is not 'as good as new'. Under these 
assumptions, a replacement policy N based on the number of repair of component 1 is considered. He determined 
the optimal replacement policy N* such that the long-run expected reward per unit time s maximized. Many 
research works have been done by Zhang [23], Stadje and Zuckerman [18], Finkelstein [6], Wang and Zhang 
[19], Lam and Zhang [10, 11]. 

Later, Wang and Zhang [20] studied a shock model with two type of failures by assuming that two kind 
of shocks in a sequence of random shocks will make the system failed, one based on the interarrival time 
between two successive shocks less than a given positive value 5 and the other based on the shock magnitude of 
single shock more than a given positive value y and the successive repair times form an increasing geometric 
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process. Under these assumptions a replacement policy N based on the number of failures of the system is 
determined such that the long-run average cost rate is minimized. Braun et.al [3] examined the increasing 
geometric process grows at most logarithmically in time, while the decreasing geometric process is almost 
certain to have a time of explosion. The a-series process grows either as a polynomial or exponential in time. It 
also noted that the geometric process doesn't satisfy a central limit theorem, while the a-series process does. 

Based on this understanding, this paper generalize Wang and Zhang [20] shock model for a repairable 
system with two-type failures model by the a -series process model. It is assumed that two kinds of shock in a 
sequence of random shocks will make the system failed, one based on the inter arrival time between two 
consecutive shocks less than a given positive value 5 (i.e. so-called '5-shock model') and the other based on the 
shock magnitude of single shock more than a given positive value 5 (i.e. so-called '5- shock model'). Assume 
further that the system after repair is 'as good as new' while the consecutive up times of the system will become 
longer and longer. Under these assumptions, by using the a-series process repair model, a replacement policy N 
based on the number of failure of the system is studied. The problem is to obtain an optimal replacement policy 
N* such that the long-run average cost per unit time is minimized. An explicit expression of the long-run average 
cost per unit time is derived, and hence an optimal policy N* can be determined analytically or numerically. 
Finally, a numerical example is given. 

For understanding the model, we considered the definition (see Braun [3]) of a- series process as follows. 
Definition 1: 

Given two random variables X and Y, if P(X>t) > P(Y>t) for all real t, then X is called stochastically 
larger than Y or Y is stochastically less than X. This is denoted by X >st Y or Y <st X respectively. 

Definition 2: 

Assume that {X k , k=l,2,. . ..}, is a sequence of independent non -negative random variables. If the distribution 
function of X k is F k (t) = F(k a t) for some a > and all k=l, 2, 3. . . then {X k , k=l, 2...} is called an a-series 
process, a is called exponent of the process. 

Obviously: 

if a>0, then {X k , k=l,2,. . ..} is stochastically decreasing, i.e, X k >st X k+1 , k=l,2,...; 
if a<0, then {X k , n=l,2,....} is stochastically increasing, i.e., X k <st X k+1 , k=l,2,; 
if a=0, then the a series process becomes a renewal process. 

II. THE A-SHOCK MODEL 

In most of the situations where random shock occurs, it can be observed that the system reliability is 
decreasing. In ordered to raise the system reliability, we consider two kinds of random shock such that the 
system is failed. Let T k and D k denote respectively the time interval between (k-l )th shock and k* shock, and the 
magnitude of k th shock. 

Let 5 and y be two positive in the system. The system will fail as soon as T k < 8 appears, k=l, 2, . . . 
and the system failure caused by T k < 8 is called the first type failure. In fact, it is a _-shock model considered by 
Wang and Zhang (2001). Similarly, the system will fail as soon as D k >y, k=l, 2, . . . appears and the system 
failure caused by D k >y is called the second type failure. In fact, it is a shock model considered by Shanthikumar 
and Sumita (1983). Remember that the system will fail as soon as both T k <8 and D k >y appears k=l, 2, . . . . 
Thus, we may as well call the system failure the third type failure. In other words, two-type failures will cause 
three cases of the system failure. It is a more extensive shock model in which we not only consider the 
magnitude of single shock but also the time interval of consecutive two shocks to the effect of the system. 

Let Si, S2, ... be shock arrival time when the system is working, and assume that 
{T k , k=l, 2, . . .} and {D k , k=l, 2, . . . } are respectively a sequence of i.i.d. non-negative random variables. 
Clearly, 

S k =T 1+ T 2 + + T k , k = l,2, (2.1) 

Let X be the working time before the system first fails, then 

Z=T 1+ T 2 + +T M . (2.2) 

Where, M is shock number of times at the system failure, and it is a random variable. Now, we give out 
some results of the system under shock model with two-type failure. 

Theorem 1: Assume that {(T k , DiJ, k = 1, 2, . . .} forms a correlated renewal sequence pair with joint 
distribution function 

F{x,y)=P{T k <x, D k <y), 
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and H(t) is the distribution function of X, then the Laplace transform of H(t) is given by 

(2.3) 

GO CO 

where F T (s) = \e~ sx F T (x) dx, F* (s, y) — \e dF T (x, y), 8 and y are 



two given positive parameters 
[See G. J. Wang & Y. L. Zhang (2005)] 

Thus, we can get the system reliability R(t) = H (t) . 
From Theorem 1 we have 

Poo = P(T>5, D<y), (2.4) 
P 10 = P(T<5, T><y), (2.5) 
Poi = P(T>5, D>y), (2.6) 
and 

P„=p(T<5, D>y), (2.7) 
denote respectively the probabilities that neither the first type failure nor the second type failure happens, only 
the first type failure happens, only the second type failure happens and both the two-type failures happen in the 
system. Clearly, p o+Pio+Poi+Pn=l- 

Theorem 2: Under the condition of Theorem 1 given by G. J. Wang & Y. L. Zhang [20] 
, the mean working time before the system fails is given by: 

FT 

EX =E[T, +T 2 + +T M ] =- . 

1 Poo (2.8) 
We make the following assumptions about the repairable system for random shock with two-type failures. 

1 . In the initial stage, the system is new. The system will be repaired as soon as it failed. The system will 
be replaced some time by a new and identical one, and the replacement time is negligible. 

2. Assume that the system after repair is 'as good as new' while the consecutive repair times of the 
system form a stochastic increasing a- series process. The time interval between the completion of the 
(k-l)-th repair and the completion of the k-th repair on the system is called the k-th cycle of the system, 
k=l, 2, . . . . 

3. Let X k l and Y k be respectively the working time and the repair time of the system incurred by i-th 

type failure in the k-th cycle, i=l, 2, 3; k=l, 2, . . . . and both the times form a decreasing a- series 
process and increasing a- series process respectively. 

4. Let {Xf , k=l, 2, . . . } is a sequence of the i.i.d. non-negative random variables distributed 

according to an exponential failure law with distribution function F(k a t) , while the distribution of 

Yjf* is assumed to be G(k^ t) , and the assume that ETj=X and E Y^ =u, is where t>0, 0<p<l, A>0, pi 
>0; i=l, 2, 3; k=l, 2, . . . . 

5. Assume that X <° , Y k (i) , i=l, 2, 3; k=l, 2, . . . are independent. 

6. Let ET = A, be the average working time. 

7. The replacement policy N based on the failure number of the system is used. 

8. Assume that the repair cost incurred by i-th type failure of the system per unit time is C r , i=l, 2, 3, 
the working reward of the system per unit time is c w , and the replacement cost of the system is C. 

III. LONG-RUN AVERAGE COST RATE UNDER POLICY N 

Now, our problem is to determine an optimal replacement policy N such that the long-run average cost 
per unit time is minimized. 

Let T] be the first replacement time, and let x k (k>l) be the time between the (k-l)-th replacement and 
the k-th replacement. Obviously, {ii, 12,13, . . .} forms a renewal process. Let C(N) be the long-run average cost 
per unit time under the replacement policy N. Because { X\, T2, T3, . . .} is a renewal process, the time interval 
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between two consecutive replacements is a renewal cycle. Then according to the renewal reward theorem (see, 
for example, Ross [13-15]), we have 

the exp ected cos? incurred in a cycle 

C(N) = . 

the exp ected length of a cycle 

Where the length of the renewal cycle is given by 



(3.1) 



w = x l + [y^ +Y l (2) / : (2) +y< 3) 7 : (3) ]+ x 2 + [y 2 (1) /<» +y 2 <2) /f ) +y 2 (3) + x N _, 



, [vffl 7(1) , V(2) 7(2) , y(3) r(3) 1 , Y 
+ l*N-l 1 AM +I N-X 1 AM + 1 AM 7 AT-1 J + A / 



(3.2) 



To evaluate EW, we only need to calculate the following results by using the property of the conditional 
expectation 

E(Y" /«) = e[ E (Y" /« I /«)] = P(/« = 1) 

E(Y™ 7f>) = £[£(7< 2 ' 7f I /<»)] = P(7f = 1) E(Y^) 

E(Y™ /?>) = E[E(7f> 7f > I 7f )] = P(7f = 1) E(y«) 
Let qj, i=l, 2, 3 be the probability of the i th type failure when failure happen, then 



^0 = ^ = p(/< 2) =i) 



E(I™) = q 3 =j P(/f>=l) - 
Clearly, q ; + q 2 +q 3 =1 ; and 
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(3.11) 

According to the assumption (6) of the model and using equations (3.9) to (3.11), the expected length 
of a renewal cycle is 
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According to equations (3.1) and (3.2), we have 
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According to equations (3.3) and (3.13), we have 
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Where 

In the next section, numerical results are provided to highlight the obtained theoretical results. 

IV. NUMERICAL RESULTS AND CONCLUSIONS: 

For the given hypothetical values of the parameters pi, P2, P3, K C r , C, (2) , C, (3) , qi, q 2 , q3 , Pi, p 2 , P 3 , C 
,and C w we determine the long run average cost per unit time from the equation (9) as follows: 
Where p 1= 20, p 2 =30, p 3 =50, ; 1=80, C r (1) =3, C r (2) =5, C, (3) =10, q,=0.45. q 2 =0.4. q 3 =0.15 , p,=-0.65. p 2 =-0.75, p 3 =- 
0.85,C=7000 ,and C w =50. 
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Table 4.1: Values of the Long-Run Average Cost Rate Under Policy N. 

V. CONCLUSIONS 

1. From the table 4.1 and graph 4.1, it is examined that the long-run average cost per unit time is minimum 
when the number of failure of the system reaches 6. i.e., C(6)=-27.751 1 at P^-0.65, p 2 =-0.75 , p 3 =- 
0.85. Thus the system should be replaced at the time of 6 th failure. 

2. The decreasing CC -series process may be more appropriate for modeling system working times while the 
increasing geometric process is more suitable for modeling repair times of the system. It assumed that the 
system after repair is not 'as good as new' and also the successive working times form a decreasing CC - 
series process, the successive repair time's form an increasing a- series process and both the processes are 
exposing to exponential failure law. Under these assumptions we study an optimal replacement policy N in 
which we replace the system when the number of failures reaches N. We determine an optimal repair 
replacement policy N* such that the long run average cost per unit time is minimized. 

3. If the repair man is familiar with the repair the successive repair times stochastically non -increasing and the 
successive working times are stochastically non -decreasing. Thus this model can also be applied for an 
improved system exposing to Weibull failure law. 
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